Meiotic drive against chromosome fusions in butterfly hybrids.

Species frequently differ in the number and structure of chromosomes they harbor, but individuals that are heterozygous for chromosomal rearrangements may suffer from reduced fitness. Chromosomal rearrangements like fissions and fusions can hence serve as a mechanism for speciation between incipient lineages, but their evolution poses a paradox. How can rearrangements get fixed between populations if heterozygotes have reduced fitness? One solution is that this process predominantly occurs in small and isolated populations, where genetic drift can override natural selection. However, fixation is also more likely if a novel rearrangement is favored by a transmission bias, such as meiotic drive. Here, we investigate chromosomal transmission distortion in hybrids between two wood white (Leptidea sinapis) butterfly populations with extensive karyotype differences. Using data from two different crossing experiments, we uncover that there is a transmission bias favoring the ancestral chromosomal state for derived fusions, a result that shows that chromosome fusions actually can fix in populations despite being counteracted by meiotic drive. This means that meiotic drive not only can promote runaway chromosome number evolution and speciation, but also that it can be a conservative force acting against karyotypic change and the evolution of reproductive isolation. Based on our results, we suggest a mechanistic model for why chromosome fusion mutations may be opposed by meiotic drive and discuss factors contributing to karyotype evolution in Lepidoptera.

A Recql5 mutant facilitates complex CRISPR/Cas9-mediated-chromosomal engineering in mouse zygotes.

Complex chromosomal rearrangements (CCRs) are often observed in clinical samples from patients with cancer and congenital diseases but are difficult to induce experimentally. Here, we report the first success in establishing animal models for CCRs. Mutation in Recql5, a crucial member of the DNA helicase RecQ family involved in DNA replication, transcription, and repair, enabled CRISPR/Cas9-mediated CCRs, establishing a mouse model containing triple fusion genes and megabase-sized inversions. Some of these structural features of individual chromosomal rearrangements use template switching and microhomology-mediated break-induced replication mechanisms and are reminiscent of the newly described phenomenon "chromoanasynthesis." These data show that Recql5-mutant mice could be a powerful tool to analyze the pathogenesis of CCRs (particularly chromoanasynthesis) whose underlying mechanisms are poorly understood. The Recql5 mutants generated in this study are to be deposited at key animal research facilities, thereby making them accessible for future research on CCRs.

Detection of 4p16.3 deletion and 11p15.5p15.4 gain in a boy by comparative genomic hybridization array: A case report.

BACKGROUND: Nonallelic homologous recombination (NAHR) of segmental duplications or low copy repeats (LCRs) result in DNA gain/loss and play an important role in the origin of genomic disorders. CASE SUMMARY: A 3-year- old boy was referred for genetic analysis. Comparative genomic hybridization array analysis revealed a loss of 3776 kb in the 4p16.3 chromosomal region and a gain of 3201 kb in the 11p15.5p15.4 chromosomal region. CONCLUSION: Genomic imbalances caused by NAHR in LCRs result in deletion and duplication syndromes.

Central nervous system embryonal tumors with EWSR1-PLAGL1 rearrangements reclassified as INI-1 deficient tumors at relapse.

PURPOSE: Central nervous system (CNS) embryonal tumors are a diverse group of malignant tumors typically affecting pediatric patients that recently have been better defined, and this paper describes evolution of a unique type of embryonal tumor at relapse. METHODS: Two pediatric patients with CNS embryonal tumors with EWSR1-PLAGL1 rearrangements treated at Arkansas Children's Hospital with histopathologic and molecular data are described. RESULTS: These two patients at diagnosis were classified as CNS embryonal tumors with EWSR1-PLAGL1 rearrangements based on histologic appearance and molecular data. At relapse both patient's disease was reclassified as atypical teratoid rhabdoid tumor (ATRT) based on loss of INI-1, presence of SMARCB1 alterations, and methylation profiling results. CONCLUSION: CNS embryonal tumors with EWSR1-PLAGL1 rearrangements acquire or include a population of cells with SMARCB1 alterations that are the component that predominate at relapse, suggesting treatment aimed at this disease component at diagnosis should be considered.

Comparative karyotype analysis of the red brocket deer (M. americana sensu lato and M. rufa) complex: evidence of drastic chromosomal evolution and implications on speciation process.

Chromosomal rearrangements are often associated with playing a role in the speciation process. However, the underlying mechanism that favors the genetic isolation associated with chromosomal changes remains elusive. In this sense, the genus Mazama is recognized by its high level of karyotype diversity among species with similar morphology. A cryptic species complex has been identified within the genus, with the red brocket deer (Mazama americana and Mazama rufa) being the most impressive example. The chromosome variation was clustered in cytotypes with diploid numbers ranging from 42 to 53 and was correlated with geographical location. We conducted an analysis of chromosome evolution of the red brocket deer complex using comparative chromosome painting and Bacterial Artificial Chromosome (BAC) clones among different cytotypes. The aim was to deepen our understanding of the karyotypic relationships within the red brocket, thereby elucidating the significant chromosome variation among closely related species. This underscores the significance of chromosome changes as a key evolutionary process shaping their genomes. The results revealed the presence of three distinct cytogenetic lineages characterized by significant karyotypic divergence, suggesting the existence of efficient post-zygotic barriers. Tandem fusions constitute the main mechanism driving karyotype evolution, following a few centric fusions, inversion X-autosomal fusions. The BAC mapping has improved our comprehension of the karyotypic relationships within the red brocket deer complex, prompting questions regarding the role of these changes in the speciation process. We propose the red brocket as a model group to investigate how chromosomal changes contribute to isolation and explore the implications of these changes in taxonomy and conservation.

Big jaw, small genome: first description of the mitochondrial genome of Odontomachus (Formicidae, Ponerinae): evolutionary implications for Ponerinae ants.

Mitochondrial DNA is a valuable tool for population genetics and evolutionary studies in a wide range of organisms. With advancements in sequencing techniques, it's now possible to gain deeper insights into this molecule. By understanding how many genes there are, how they're organized within the molecule, identifying the presence of spacers, and analyzing the composition of the D-Loop, we can better grasp the rearrangements that play a crucial role in the evolutionary dynamics of mitochondrial DNA. Additionally, phylogenetic analyses benefit significantly from having access to a larger pool of mtDNA genes. This wealth of genetic information allows for the establishment of evolutionary relationships with greater accuracy than ever before, providing a more robust framework than analyses based on a limited number of genes. Studies on mitogenomes belonging to the family Formicidae have proven promising, enabling the identification of gene rearrangements and enhancing our understanding of the internal relationships within the group. Despite this, the number of mitogenomes available for the subfamily Ponerinae is still limited, and here we present for the first time the complete mitogenome of Odontomachus. Our data reveal a gene duplication event in Formicidae, the first involving trnV, and new gene arrangements involving the trnM-trnI-trnQ and trnW-trnC-trnY clusters, suggesting a possible synapomorphy for the genus. Our phylogenetic analysis using the PCGs available for Formicidae supports the monophyly of the subfamily Ponerinae and sheds light on the relationship between Odontomachus and Pachycondyla.

How structural variants shape avian phenotypes: Lessons from model systems.

Despite receiving significant recent attention, the relevance of structural variation (SV) in driving phenotypic diversity remains understudied, although recent advances in long-read sequencing, bioinformatics and pangenomic approaches have enhanced SV detection. We review the role of SVs in shaping phenotypes in avian model systems, and identify some general patterns in SV type, length and their associated traits. We found that most of the avian SVs so far identified are short indels in chickens, which are frequently associated with changes in body weight and plumage colouration. Overall, we found that relatively short SVs are more frequently detected, likely due to a combination of their prevalence compared to large SVs, and a detection bias, stemming primarily from the widespread use of short-read sequencing and associated analytical methods. SVs most commonly involve non-coding regions, especially introns, and when patterns of inheritance were reported, SVs associated primarily with dominant discrete traits. We summarise several examples of phenotypic convergence across different species, mediated by different SVs in the same or different genes and different types of changes in the same gene that can lead to various phenotypes. Complex rearrangements and supergenes, which can simultaneously affect and link several genes, tend to have pleiotropic phenotypic effects. Additionally, SVs commonly co-occur with single-nucleotide polymorphisms, highlighting the need to consider all types of genetic changes to understand the basis of phenotypic traits. We end by summarising expectations for when long-read technologies become commonly implemented in non-model birds, likely leading to an increase in SV discovery and characterisation. The growing interest in this subject suggests an increase in our understanding of the phenotypic effects of SVs in upcoming years.

An Optimized Genotyping Workflow for Identifying Highly SCRaMbLEd Synthetic Yeasts.

Synthetic Sc2.0 yeast strains contain hundreds to thousands of loxPsym recombination sites that allow restructuring of the Saccharomyces cerevisiae genome by SCRaMbLE. Thus, a highly diverse yeast population can arise from a single genotype. The selection of genetically diverse candidates with rearranged synthetic chromosomes for downstream analysis requires an efficient and straightforward workflow. Here we present loxTags, a set of qPCR primers for genotyping across loxPsym sites to detect not only deletions but also inversions and translocations after SCRaMbLE. To cope with the large number of amplicons, we generated qTagGer, a qPCR genotyping primer prediction tool. Using loxTag-based genotyping and long-read sequencing, we show that light-inducible Cre recombinase L-SCRaMbLE can efficiently generate diverse recombination events when applied to Sc2.0 strains containing a linear or a circular version of synthetic chromosome III.

Clinical outcomes in carriers of insertional translocation: a retrospective analysis of comprehensive chromosome screening results.

Objective: To evaluate the clinical outcomes in the carriers of insertional translocation (IT). Design: Retrospective case series. Setting: University-affiliated reproductive medical center. Patients: Twenty-three couples with ITs. Intervention: No direct interventions were involved; however, this study included patients who underwent preimplantation genetic testing for structural chromosomal rearrangements (PGT-SR). Main Outcome Measure: Outcome of preimplantation genetic testing for structural chromosomal rearrangements and percentage of blastocysts available for transfer. Results: Among 23 IT carriers, 15 were simple interchromosome ITs (type A), 3 were intrachromosome IT carriers (type B), and 5 were interchromosome IT carriers combined with other translocations (type C). A total of 190 blastocysts from 30 cycles were biopsied, 187 embryos were tested successfully, and only 57 blastocysts (30.5%) from 21 patients were available for transfer (normal or balanced). The unbalanced rearrangement rate of type C was 79.2% (42/53), and the proportion of type A was 50.0% (57/114), which was significantly higher than that of type B (5%, 1/20). In type A, the probability of embryos harboring unbalanced rearrangement in female carriers was 56.0% (51/91), which was higher than that in male carriers (26.1%, 6/23). Furthermore, the haploid autosomal length value of the inserted fragment was correlated linearly with the incidence of abnormal embryos. In type A gametes, most gametes produced by 2:2 separation without crossover, and no 3:1 separation gamete was observed. Conclusions: The chance of identifying normal or balanced and mosaic blastocysts per mature oocytes in patients with ITs are 16.6% (67/404). Greater IT complexity results in fewer transferable embryos. For simple interchromosome ITs, female carriers and those with higher haploid autosomal length values have a higher risk of producing embryos with unbalanced rearrangement.

Optical Genome Mapping as a Potential Routine Clinical Diagnostic Method.

Chromosome analysis (CA) and chromosomal microarray analysis (CMA) have been successfully used to diagnose genetic disorders. However, many conditions remain undiagnosed due to limitations in resolution (CA) and detection of only unbalanced events (CMA). Optical genome mapping (OGM) has the potential to address these limitations by capturing both structural variants (SVs) resulting in copy number changes and balanced rearrangements with high resolution. In this study, we investigated OGM's concordance using 87 SVs previously identified by CA, CMA, or Southern blot. Overall, OGM was 98% concordant with only three discordant cases: (1) uncalled translocation with one breakpoint in a centromere; (2) uncalled duplication with breakpoints in the pseudoautosomal region 1; and (3) uncalled mosaic triplication originating from a marker chromosome. OGM provided diagnosis for three previously unsolved cases: (1) disruption of the SON gene due to a balanced reciprocal translocation; (2) disruption of the NBEA gene due to an inverted insertion; (3) disruption of the TSC2 gene due to a mosaic deletion. We show that OGM is a valid method for the detection of many types of SVs in a single assay and is highly concordant with legacy cytogenomic methods; however, it has limited SV detection capabilities in centromeric and pseudoautosomal regions.

Reinstatement of the fission yeast species Schizosaccharomyces versatilis Wickerham et Duprat, a sibling species of Schizosaccharomyces japonicus.

Schizosaccharomyces japonicus Yukawa et Maki (1931) and Schizosaccharomyces versatilis Wickerham et Duprat (1945) have been treated as varieties of S. japonicus or as conspecific, based on various approaches including mating trials and nDNA/nDNA optical reassociation studies. However, the type strains of S. japonicus and S. versatilis differ by five substitutions (99.15% identity) and one 1-bp indel in the sequences of the D1/D2 domain of the 26S rRNA gene, and 23 substitutions (96.3% identity) and 31-bp indels in the sequences of internal transcribed spacer (ITS) of rRNA, suggesting that they may not be conspecific. To reassess their taxonomic status, we conducted mating trials and whole-genome analyses. Mating trials using the type strains showed a strong but incomplete prezygotic sterility barrier, yielding interspecies mating products at two orders of magnitude lower efficiency than intraspecies matings. These mating products, which were exclusively allodiploid hybrids, were unable to undergo the haplontic life cycle of the parents. We generated chromosome-level gap-less genome assemblies for both type strains. Whole genome sequences yielded an average nucleotide identity (ANI) of 86.4%, indicating clear separation of S. japonicus and S. versatilis. Based on these findings, we propose the reinstatement of S. versatilis as a distinct species (holotype strain: CBS 103T and ex-types: NRRL Y-1026, NBRC 1607, ATCC 9987, PYCC 7100; Mycobank no.: 847838).

Gene Rearrangements in the Mitochondrial Genome of Gonatopsis borealis and Onychoteuthis compacta Reveal Their Phylogenetic Implications for Oegopsida.

The complete mitochondrial genome provides crucial information for comprehending gene rearrangement, molecular evolution, and phylogenetic analysis. Here, we have determined the complete mitogenome sequence of Gonatopsis borealis and Onychoteuthis compacta for the first time. Their genome sizes were 20,148 bp and 20,491 bp, respectively, including 18 protein-coding genes, COI-COIII, ATP6, and ATP8 are duplicated, 23 transfer RNA genes, and 2 ribosomal RNA (rRNA) genes (12S and 16S rRNA). Specifically, the overall A+T content is 70.69% and 72.67%. It shows a significant AT bias. The whole mitogenomes indicate positive AT skew (0.070 and 0.062). Furthermore, the gene order has been rearranged within Oegopsida. The tandem duplication random loss model was determined as most likely to explain the observed gene rearrangements. Phylogenetic analysis was performed, and the result tree was found to be consistent with the morphological identification classification. Estimation of divergence time for 35 species showed that the main differentiation of Oegopsida occurred in 140.70 Mya. These results will help to better understand the gene rearrangements and evolution of G. borealis and O. compacta and lay a foundation for further phylogeny genetic studies of Oegopsida.

NG2 Molecule Expression in Acute Lymphoblastic Leukemia B Cells: A Flow-Cytometric Marker for the Rapid Identification of KMT2A Gene Rearrangements.

Background: B-lineage acute lymphoblastic leukemias (B-ALL) harboring rearrangements of the histone lysine [K]-Methyltransferase 2A (KMT2A) gene on chromosome 11q23 (KMT2A-r) represent a category with dismal prognosis. The prompt identification of these cases represents an urgent clinical need. Considering the correlation between rat neuron glial-antigen 2 (NG2) chondroitin-sulfate-proteoglycan molecule expression and KMT2A-r, we aimed to identify an optimized cytofluorimetric diagnostic panel to predict the presence of KMT2A-r. Materials and Methods: We evaluated 88 NG2+ B-ALL cases identified with an NG2 positivity threshold >10% from a cohort of 1382 newly diagnosed B-ALLs referred to the Division of Hematology of 'Sapienza' University of Rome. Results: Eighty-five of 88 (96.6%) NG2+ B-ALLs harbored KMT2A-r and were mainly pro-B ALL (77/85; 91%). Only 2 B-ALLs with KMT2A-r showed NG2 expression below 10%, probably due to the steroid therapy administered prior to cytofluorimetric analysis.Compared to KMT2A-r-cases, KMT2A r+ B-ALLs showed a higher blast percentage, significantly higher mean fluorescence intensity (MFI) of CD45, CD38, and CD58, and significantly lower MFI of CD34, CD22, TdT, and CD123.The study confirmed differences in CD45, CD34, CD22, and TdT MFI within the same immunologic EGIL group (European Group for the immunological classification of leukemias), indicating no influence of the B-ALLs EGIL subtype on the KMT2A-r+ B-ALLs immunophenotype. Conclusions: Our data demonstrate the association between NG2 and KMT2A-r in B-ALLs identify a distinctive immunophenotypic pattern, useful for rapid identification in diagnostic routines of these subtypes of B-ALLs with a poor prognosis that benefits from a specific therapeutic approach.

Nanopore sequencing with T2T-CHM13 for accurate detection and preventing the transmission of structural rearrangements in highly repetitive heterochromatin regions in human embryos.

BACKGROUND: Structural rearrangements in highly repetitive heterochromatin regions can result in miscarriage or foetal malformations; however, detecting and preventing the transmission of these rearrangements has been challenging. Recently, the completion of sequencing of the complete human genome (T2T-CHM13) has made it possible to accurately characterise structural rearrangements in these regions. We developed a method based on T2T-CHM13 and nanopore sequencing to detect and block structural rearrangements in highly repetitive heterochromatin sequences. METHODS: T2T-CHM13-based "Mapping Allele with Resolved Carrier Status" was performed for couples who carry structural rearrangements in heterochromatin regions. Using nanopore sequencing and the T2T-CHM13 reference genome, the precise breakpoints of inversions and translocations close to the centromere were detected and haplotypes were constructed using flanking single-nucleotide polymorphisms (SNPs). Haplotype linkage analysis was then performed by comparing consistent parental SNPs with embryonic SNPs to determine whether the embryos carried hereditary inversions or balanced translocations. Based on copy number variation and haplotype linkage analysis, we transplanted normal embryos, which were further verified by an amniotic fluid test. RESULTS: To validate this approach, we used nanopore sequencing of families with inversions and reciprocal translocations close to the centromere. Using the T2T-CHM13 reference genome, we accurately detected inversions and translocations in centromeres, constructed haplotypes and prevented the transmission of structural rearrangements in the offspring. CONCLUSIONS: This study represents the first successful application of T2T-CHM13 in human reproduction and provides a feasible protocol for detecting and preventing the transmission of structural rearrangements of heterochromatin in embryos.

The Landscape and Predicted Roles of Structural Variants in Fusarium Graminearum Genomes.

Structural rearrangements, such as inversions, translocations, duplications, and large insertions and deletions, are large-scale genomic variants that can play an important role in shaping phenotypic variation and in genome adaptation and evolution. We used chromosomal-level assemblies from eight Fusarium graminearum isolates to study structural variants and their role in fungal evolution. We generated the assemblies of four of these genomes after Oxford Nanopore sequencing. A total of 87 inversions, 159 translocations, 245 duplications, 58,489 insertions and 34,102 deletions were detected. Regions of high recombination rate are associated with structural rearrangements, and a significant proportion of inversions, translocations, and duplications overlap with the repeat content of the genome, suggesting recombination and repeat elements are major factors in the origin of structural rearrangements in F. graminearum. Large insertions and deletions introduce presence-absence polymorphisms for many genes, including secondary metabolite biosynthesis cluster genes and predicted effectors genes. Translocation events were found to be shuffling predicted effector-rich regions of the genomes and are likely contributing to the gain and loss of effectors facilitated by recombination. Breakpoints of some structural rearrangements fall within coding sequences and are likely altering the protein products. Structural rearrangements in F. graminearum thus have an important role to play in shaping pathogen-host interactions and broader evolution through genome reorganization, the introduction of presence-absence polymorphisms, and changing protein products and gene regulation.

Fifth Edition of the World Health Organization Classification of Tumors of the Hematopoietic and Lymphoid Tissues: Acute Lymphoblastic Leukemias, Mixed-Phenotype Acute Leukemias, Myeloid/Lymphoid Neoplasms With Eosinophilia, Dendritic/Histiocytic Neoplasms, and Genetic Tumor Syndromes.

This manuscript represents a review of lymphoblastic leukemia/lymphoma (acute lymphoblastic leukemia/lymphoblastic lymphoma), acute leukemias of ambiguous lineage, mixed-phenotype acute leukemias, myeloid/lymphoid neoplasms with eosinophilia and defining gene rearrangements, histiocytic and dendritic neoplasms, and genetic tumor syndromes of the 5th edition of the World Health Organization Classification of Tumors of the Hematopoietic and Lymphoid Tissues. The diagnostic, clinicopathologic, cytogenetic, and molecular genetic features are discussed. The differences in comparison to the 4th revised edition of the World Health Organization classification of hematolymphoid neoplasms are highlighted.

Structural Changes in the Carbon Sphere of a Dirhodium Complex Induced by Redox or Deprotonation Reactions.

A carbon-rich molecule is synthesized, which mainly contains conjugated sp2 and sp hybridized carbon centers. Alkenyl and alkynyl binding sites are arranged such that this compound serves as ligand to a binuclear metal unit with a RhI─RhI bond. Furthermore, CH units are placed in proximity to the metal centers. The dicationic complex [Rh2(bipy)2{Ph2Ptrop(C≡CCy)2}]2+(OTf-)2 allows to study possible responses of the carbon-framework to redox reactions as well as deprotonation reactions. All products are, whenever possible, characterized by X-ray diffraction (XRD) methods, NMR and EPR spectroscopy as well as electrochemical methods. It is shown that the carbon skeleton of the ligand framework undergoes C─C bond rearrangement reactions of remarkable diversity. In combination with DFT (density functional theory) studies, these results allow to gain insight into the electronic structure changes caused by metal sites in a carbon-rich environment, which may be of relevance for the properties of metal particles on carbon support materials when they are exposed to hydrogen, electrons, or protons.

Uncoupling programmed DNA cleavage and repair scrambles the Paramecium somatic genome.

In the ciliate Paramecium, precise excision of numerous internal eliminated sequences (IESs) from the somatic genome is essential at each sexual cycle. DNA double-strands breaks (DSBs) introduced by the PiggyMac endonuclease are repaired in a highly concerted manner by the non-homologous end joining (NHEJ) pathway, illustrated by complete inhibition of DNA cleavage when Ku70/80 proteins are missing. We show that expression of a DNA-binding-deficient Ku70 mutant (Ku70-6E) permits DNA cleavage but leads to the accumulation of unrepaired DSBs. We uncoupled DNA cleavage and repair by co-expressing wild-type and mutant Ku70. High-throughput sequencing of the developing macronucleus genome in these conditions identifies the presence of extremities healed by de novo telomere addition and numerous translocations between IES-flanking sequences. Coupling the two steps of IES excision ensures that both extremities are held together throughout the process, suggesting that DSB repair proteins are essential for assembly of a synaptic precleavage complex.

Uncovering rearrangements in the Tibetan antelope via population-derived genome refinement and comparative analysis with homologous species.

Introduction: The Tibetan antelope (Pantholops hodgsonii) is a remarkable mammal thriving in the extreme Qinghai-Tibet Plateau conditions. Despite the availability of its genome sequence, limitations in the scaffold-level assembly have hindered a comprehensive understanding of its genomics. Moreover, comparative analyses with other Bovidae species are lacking, along with insights into genome rearrangements in the Tibetan antelope. Methods: Addressing these gaps, we present a multifaceted approach by refining the Tibetan Antelope genome through linkage disequilibrium analysis with data from 15 newly sequenced samples. Results: The scaffold N50 of the refined reference is 3.2 Mbp, surpassing the previous version by 1.15-fold. Our annotation analysis resulted in 50,750 genes, encompassing 29,324 novel genes not previously study. Comparative analyses reveal 182 unique rearrangements within the scaffolds, contributing to our understanding of evolutionary dynamics and species-specific adaptations. Furthermore, by conducting detailed genomic comparisons and reconstructing rearrangements, we have successfully pioneered the reconstruction of the X-chromosome in the Tibetan antelope. Discussion: This effort enhances our comprehension of the genomic landscape of this species.

The complete mitochondrial genome of Penicillium expansum: Insights into the fungal evolution and phylogeny.

Penicillium expansum is an important phytopathogenic fungus that causes blue mold disease. In this study, the novel mitochondrial genome of P. expansum was sequenced, assembled, annotated, and compared with the previously published Penicillium mitogenomes. P. expansum mitogenome is composed of circular DNA molecules with a genome size of 25,496 bp. It encodes 16 protein-encoding genes (PCGs), two rRNA genes, and 25 tRNA genes. Comparative analysis with six other Penicillium species revealed that gene length, GC content, AT skew, and GC skew were variable among the core protein-coding genes. The Penicillium species' gene synteny analysis identified several gene rearrangements. Among the core 15 PCGs, atp8 had the lowest K2P genetic distance, which shows that this gene is highly conserved. The Ka/Ks value of most PCGs was less than 1, which shows that these genes have undergone purifying selection. Phylogenetic analysis based on 14 concatenated core mitochondrial genes revealed that P. expansum shares a close relationship with P. solitum. This study served as a first report on the complete mitochondrial genome of P. expansum and its comparative analysis that will contribute to population genetics and rapid evolutionary studies among Penicillium species.